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This  is  a  report  of  progress  during  the  past  year.  The  work  was  performed  in 
three  areas  with  a  long  term  goal  understanding  the  formation  and  maintenance 
of  electrostatic  fields  in  the  earth's  magnetosphere.  The  entry  of  low  energy 
charged  particles  into  a  magnetically  closed  magnetosphere  has  been  examined 
in  some  detail.  Entry  is  permitted  because  of  the  non-uniform  nature  of  the 
magnetic  field  over  the  magnetopause  surface.  Electrostatic  fields  may  be 
formed  across  the  tail  of  the  magnetosphere  because  of  the  different  "entry 
efficiencies"  of  protons  and  electrons.  The  consequences  of  this  nartirlex 
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entry  mechanism  for  the  plasma  sheet,  plasma  mantle,  and  boundary  plasmas 
in  the  magnetosphere  are  examined.  The  mathematics  of  particle  entry 
was  investigated  in  a  one-dimensional  boundary  using  both  kinetic  theory 
bulk  MHD  parameters.  From  our  participation  in  the  6th  Coordinated 
Data  Analysis  Workshop,  we  have  determined  that  at  least  during  disturbed 
magnetic  conditions,.  Currents  persist  near  geoshychronous  orbit  in  the 
nightime  region  which  are  presently  not  included  in  our  dynamic  magnetic 
field  models.  These  currents  are  probably  associated  with  the  field 
aligned  currents  which  close  in  the  ionosphere  near  auroral  latitudes. 
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1.0  INTRODUCTION 

The  work  of  the  McDonnell  Douglas  Group  supported  by  ONR  support  has  always 
been  focused  on  magnetospheric  electric  fields.  In  our  early  work,  we  found 
that  in  order  to  quantitatlvly  represent  induced  electric  fields,  we  had  to 
first  accurately  model  the  magnetospheric  magnetic  field.  This  took  several 
years  and  resulted  in  models  that  are  now  used  worldwide  to  represent  the 
extension  of  the  geomagnetic  field  into  space.  More  recently  we  have 
attempted  to  quantitatively  represent  the  electrostatic  fields  present  in  the 
magnetosphere.  The  charge  separation  responsible  for  this  field  is 
Inextricably  tied  to  the  microphysics  of  the  magnetopause  and  charged  particle 
entry  into  a  magnetically  closed  magnetosphere.  Thus,  again,  we  found  that  we 
must  first  understand  and  model  another  magnetospheric  feature— the  entry  and 
motion  of  low  energy  particles.  In  particular  we  have  been  Interested  in  the 
formation  of  the  cross  tall  electric  field  which  Is  of  electrostatic  origin. 

Our  interest  in  this  problem  led  us  to  reexamine  the  classical  pressure 
balance  condition  which  has  been  used  historically  to  determine  the  shape  of 
the  magnetosphere  and  the  topology  of  the  magnetospheric  magnetic  field.  The 
pressure  balance  equation  states  simply  that  the  transition  from  the  extension 
of  the  geomagnetic  field  to  the  region  of  interplanetary  space  will  be  found 
where  the  energy  density  (or  pressure)  of  the  geomagnetic  field  is  equal  to 
the  pressure  caused  by  the  streaming  solar  wind  plasma.  Implicit  In  the  use 
of  the  pressure  balance  equation  is  the  assumption  that  the  charged  particles 
incident  on  the  magnetic  field  are  reflected  specularly;  that  is,  in  a  mirror 
like  fashion  such  that  a  particle's  angle  of  incidence  is  exactly  equal  to  its 
angle  of  reflection.  This  assumption  implies  that  the  magnetic  field  In  the 
vicinity  of  the  interaction  with  the  Incident  particle  Is  uniform.  We  noted 
as  early  as  our  1974  papers  that  this  was  Indeed  not  the  case  and  that 
gradients  in  the  magnetic  field  parallel  to  the  magnetopause  might  be 
responsible  for  permitting  the  entry  of  a  fraction  of  the  Incident  particles, 
the  size  of  the  fraction  being  determined  by  the  angle  between  the  incident 
particle  direction  and  the  magnetopause  surface. 
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In  the  past  year  and  a  half,  we  have  explored  this  more  realistic 
representation  of  pressure  balance  with  the  goal  of  understanding  the  entry  of 
the  magnetosheath  plasma  to  the  magnetosphere.  We  are  now  In  a  position  to 
qualitatively  explain  the  exl stance  of  the  plasma  sheet,  boundary  lay^r  and 
plasma  mantle  regions  of  the  magnetosphere.  We  are  also,  of  course, 
interested  in  understanding  better  the  microphysics  of  the  magnetopause  with 
this  new  pressure  balance  formalism.  We  expect  that  this  work  will  also  lead 
to  a  quantitative  understanding  of  the  formation  of  the  cross  tall  electric 
field.  (This  relates  back  to  our  long  term  goal  of  quantitatively 
representing  the  total  electromagnetic  environment  for  charged  particles 
throughout  the  earth's  magnetosphere. ) 

This  report  continues  with  3  technical  sections.  We  first  discuss 
qualitatively  the  entry  of  plasma  into  the  magnetosphere  using  a  realistic 
pressure  balance  formalism.  Although  most  of  our  results  to  date  on  plasma 
entry  are  qualitative,  we  feel  that  we  can  explain  the  presence  of  most  of  the 
magnetosphere's  plasma  regions  In  terms  of  this  entry  mechanism.  In  the 
second  section  we  discuss  work  performed  in  the  last  year  on  the  microphysics 
of  the  magnetopause  that  is  tied  to  our  understanding  of  the  entry  mechanism 
discussed  in  the  previous  section.  Although  some  interesting  work  has  been 
performed,  our  tentative  conclusion  on  the  microphysics  problem  is  that  it 
will  require  more  resources  to  solve  and  that  it  must  be  solved  in  advance  of 
our  work  on  the  representation  of  the  cross  tail  electrostatic  field.  The 
last  section  provides  up  to  date  information  on  our  participation  in  the 
Coordinated  Data  Analysis  Workshop  activities  with  emphasis  on  CDAW-6  which  is 
In  progress.  These  workshops  have  provided  us  much  valuable  quantitative 
information  on  the  accuracy  of  our  magnetic  field  models  and  on  the  dynamic 
properties  of  the  magnetosphere  generally. 
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2.0  THE  ENTRY  OF  LOW  ENERGY  PARTICLES  INTO  A  MAGNETICALLY  CLOSED  MAGNETOSPHERE 


This  section  is  divided  into  two  parts.  In  the  first  we  discuss  quantitative 
work  that  has  been  done  in  the  past  year.  In  it  we  also  discuss  problems 
regarding  the  data  base  (in  particular  limited  Information  on  particle 
distribution  functions  In  the  magnetosheath  region)  and  the  Impact  that  has  on 
future  quantitative  work.  In  the  second  section  qualitative  aspects  on 
particle  entry  into  the  magnetosphere  are  discussed  generally  in  terms  of  our 
nonspecular  reflection  model  of  the  pressure  balance  condition. 

2.1  Quantitative  Work 

We  have  explained  in  earlier  reports  our  quantitative  procedures  for 
determining  charged  particle  entry  through  the  magnetopause.  This  Is  done  by 
computer  and  involves  several  codes.  First  an  accurate  description  of  the 
magnetic  field  topology  at  the  magnetopause  is  required  for  the  quantitative 
representation  of  gradients  in  the  magnetic  field.  It  Is  the  gradient  in  the 
field  that  causes  the  particles  to  not  reflect  specularly  from  the  magnetic 
field,  allowing  some  of  them  to  remain  in  the  magnetosphere.  This  process  Is 
shown  in  Figure  1.  Another  code  calculates  the  Lorentz  force  on  the  charged 
particle  and  Is  used  to  determine  the  trajectory  of  the  particle  In  the 
magnetospheric  magnetic  field.  At  a  given  point  on  the  magnetopause, 
particles  are  allowed  to  Impact  with  varying  angles  to  the  surface  normal. 

Last  year  we  developed  a  graphic  representation  for  the  “entry  cone",  that  Is, 
the  solid  angle  region  through  which  particles  may  Impact  the  magnetopause  and 
then  remain  In  the  magnetosphere  due  to  the  non-uniform  magnetic  field.  An 
example  of  this  representation  Is  shown  in  Figure  2.  Note  that  the  "entry 
fraction"  depends  on  the  position  on  the  magnetopause,  the  particle  type 
(i.e.,  proton,  electron,  helium  nucleus),  and  the  the  particle  energy.  The 
entry  fraction  Is  the  ratio  of  the  solid  angle  for  entering  particles  to  the 
total  solid  angle  at  the  point  in  question.  It  also  represents  that  fraction 
of  Isotropic  magnetosheath  plasma  particles  that  would  enter  the 
magnetosphere.  Since  the  magnetosheath  plasma  possesses  an  anti-solar  bulk 
flow  direction,  the  actual  percent  of  Incident  plasma  particles  that  enter  the 


magnetosphere  at  any  point  will  be  somewhat  larger  than  the  entry  fraction 
numbers  discussed  in  this  section.  For  a  more  realistic  representation  of  the 
magnetosheath  plasma,  particle  distribution  functions  must  be  used. 

The  representation  of  the  solid  angle  for  entering  particles  shown  in  Figure  2 
is  of  limited  usefulness  because  it  discusses  particle  entry  only  at  a  point 
on  the  magnetopause.  This  work  was  extended  to  many  points  over  the  tail  of 

the  magnetosphere  in  order  to  provide  a  picture  of  particle  entry  as  a 

function  of  location  on  the  magnetopause  and  in  terms  of  particle  energy. 

This  discription  is  provided  in  Figures  3  through  5  which  show  particle  entry 
into  the  tail  of  the  magnetosphere  for  one,  five  and  ten  kilovolt  protons  at 
different  distances  down  the  tail.  It  is  noted  that  the  entry  fraction  is 
quite  sizable  near  the  equator  but  quickly  diminishes  at  off  equatorial 
locations.  This  is  clearly  consistent  with  observations  of  the  plasma  sheet 
which  persists  across  the  equatorial  region  of  the  tail  with  a  maximum 
thickness  of  about  8  earth  radii  near  the  flanks  (i.e.,  the  region  in  the 

equatorial  tail  near  the  magnetopause  boundary).  Note  that  for  10  kilovolt 

protons  the  entry  fraction  is  almost  .3,  that  is,  approximately  30 %  of  the 
incident  protons  at  a  distance  of  40  earth  radii  down  the  tall  enter  the 
magnetosphere.  We  note  also  that  because  of  the  magnetic  field  topology, 
particles  more  than  20  degrees  above  the  equator  in  the  tail  region,  even  if 
they  do  enter,  flow  along  magnetic  field  lines  in  the  anti-solar  direction. 
Thus  the  only  portions  of  the  tail  that  are  traversed  by  charged  particles 
originating  in  the  magnetosheath  region  are  in  the  plasma  sheet  and  in  a 
narrow  boundary  layer  that  surrounds  the  northern  and  southern  extremities  of 
the  cross  section  of  the  tail.  The  regions  in  between  (the  lobes)  are 
characterized  by  the  persistence  of  a  magnetic  field  and  a  very  low  plasma 
density.  This  is  because  the  lobe  regions  are  not  connected  to  the 
magnetosheath  particle  source  by  magnetic  field  topology.  More  is  said  on 
this  subject  in  the  section  on  qualitative  results. 

As  a  corollary  exercise,  we  have  quantitatively  examined  the  implications  of 
this  particle  entry  on  the  pressure  balance  equation,  the  magnetopause  shape, 
and  the  change  in  the  strength  of  the  magnetic  field  at  the  location  of  the 
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magnetopause.  When  25 4  of  the  plasma  incident  on  the  region  of  the 
magnetopause  enters  the  magnetosphere,  the  pressure  at  the  magnetopause  is 
over  estimated  by  12-1/24.  (The  pressure  balance  condition  assumes  a  factor 
of  2  for  the  Incident  charged  particles  since  they  are  reflected.  Entering 
particles  will  still  create  pressure  In  the  boundary  but  only  half  as  much  as 
the  reflected  particles).  This  change  in  pressure  will  act  to  Increase  the 
geocentric  distance  to  the  magnetopause  by  about  34  and  act  to  decrease  the 
magnetic  field  there  by  approximately  64.  Thus  the  changes  in  magnetopause 
shape  and  magnetic  topology  are  notlcable  but  minor  and  in  the  proper 
direction.  That  is,  they  make  the  tall  slightly  more  tall -like  (a  non  dipolar 
field  geometry)  than  shown  in  existing  magnetic  field  models.  We  therefore 
conclude  that  this  entry  mechanism  produces  only  minor  changes  and  that  it 
continues  to  be  appropriate  to  use  the  old  pressure  balance  formalism  to 
determine  magnetopause  shape  and  field  topology. 

Again,  however,  we  wish  to  emphasize  that  by  relaxing  the  specular  reflection 
condition  in  the  pressure  balance  equation,  we  are  able  to  explain  the  entry 
of  low  energy  particles  into  a  magnetically  closed  magnetosphere.  Low  energy 
particle  entry  through  the  magnetopause  occurs  virtually  at  all  locations. 
Magnetic  field  topology  within  the  magnetosphere  determines  the  subsequent 
motion  of  the  particles.  It  Is  only  in  the  equatorial  region  of  the  tail  (the 
plasma  sheet  region)  that  magnetosheath  plasma  has  direct  access  to  regions 
deep  in  the  magnetosphere.  More  is  said  concerning  the  supply  of 
magnetosheath  plasma  to  the  magnetosphere  in  the  following  section. 

2.2  The  Supply  of  Magnetosheath  Plasma  to  the  Magnetosphere 

Although  much  work  to  date  has  been  don*4  on  the  quantitative  aspects  of  low 
energy  charged  particle  entry,  problems  remain  with  the  quantitative 
calculation  of  magnetospheric  plasma  densities  based  on  magnetosheath  plasma 
values  because  of  the  small  amount  of  data  on  magnetosheath  plasma 
properties.  We  feel  however,  that  it  is  important  to  discuss  qualitatively 
those  features  of  magnetospheric  plasmas  that  can  be  explained  In  terms  of  low 
energy  particle  entry  through  the  magnetopause.  We  believe  these  Include  the 
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existence  of  all  quiet  time  plasma  features  found  In  the  mid  magnetosphere  and 
the  tail.  Near  synchronous  orbit  the  ring  current  and  the  plasmasphere  are 
populated  by  particles  originating  in  the  ionosphere  and  by  particles  that 
have  been  accelerated  inward  from  the  tail  region  during  dynamic 
magnetospheric  events.  Their  existence  thus  cannot  be  explained  directly  in 
terms  of  a  magnetosheath  source  plasma. 

Thus,  in  order  to  determine  which  regions  of  the  magnetosphere  are  populated 
by  particles  from  the  magnetosheath  it  is  necessary  to  understand  the  magnetic 
control  of  particle  motions  within  the  magnetosphere. 

To  discuss  particle  motions  in  the  magnetosphere,  it  is  appropriate  to  begin 

with  the  equatorial  region.  Particles  with  90  degree  pitch  angles  (their 

velocity  vector  perpendicular  to  the  direction  of  the  magnetic  field)  will 

remain  in  the  equatorial  plane.  Furthermore,  pa"ticles  in  a  static  magnetic 

field  configuration  will  move  along  lines  of  constant  total  magnetic  field. 

Contours  of  constant  magnetic  field  are  shown  in  Figure  6.  It  is  seen  that 

there  is  a  region  in  the  inner  magnetosphere  where  the  contours  do  not 

intersect  the  magnetopause.  It  is  this  region  and  the  lobes  of  the  tall  that 

are  not  directly  linked  to  the  magnetopause  and  therefore  cannot  directly 

receive  plasma  from  the  magnetosheath  region.  Throughout  the  remainder  of  the 

equatorial  magnetosphere,  some  of  the  particles  impacting  the  magnetopause 

will  gain  access  to  the  magnetosphere  and  drift  across  and  out  the  other 

side.  As  shown  in  Figures  3  through  5,  the  number  of  incident  particles 

entering  the  magnetosphere  increases  with  distance  down  the  tail.  Also,  since 

the  magnetic  field  is  uniform  at  the  subsolar  point  (the  intersection  of  the 

sun-earth  line  with  the  magnetopause)  it  is  expected  that  the  entry  fraction 

near  the  nose  of  the  magnetosphere  will  be  quite  small.  Thus,  in  the 

equatorial  region,  there  exists  an  inner  magnetosphere  area  where  no 

magnetosheath  particles  have  direct  access.  Moving  along  the  magnetopause 

tailward,  the  number  of  entering  particles  increases  gradually  so  that  at 

about  x  =  -8  R  the  inner  edge  of  the  plasma  sheet  is  encountered.  Note 
sm  e 

that  for  90  degree  pitch  angle  particles  the  inner  edge  of  the  plasma  sheet 
can  be  no  closer  than  about  7  Rg  from  the  center  of  the  earth  during 
undisturbed  magnetic  conditions. 
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When  the  discussion  of  particle  entry  Is  expanded  to  include  non-equatori al 
particles  (i.e.,  those  with  pitch  angles  other  than  90  degrees)  things  become 
more  complicated  since,  in  addition  to  the  gradient  drift,  the  particles  also 
move  up  and  down  along  field  lines.  However,  it  is  the  non-equatori al 
particles  in  the  magnetosphere  that  populate  the  plasma  sheet,  the  plasma 
mantle,  and  boundary  layer  plasmas.  We  believe  that  the  existence  of  these 
different  plasma  populations  in  the  magnetosphere  can  be  attributed  to  one 
phenomenon;  the  entry  of  low  energy  plasma  into  a  magnetically  closed 
magnetosphere.  Also  note  that  although  the  following  discussion  is 
qualitative  with  regard  to  charged  particle  densities,  the  comments  on  the 
locations  of  these  various  plasmas  within  the  magnetosphere  are  quantitative. 

A  cross  section  of  the  magnetospheric  tail  is  shewn  schematically  in  Figure  7 
where  several  observed  features  are  noted.  The  plasma  sheet  persists  through 
the  center  of  the  tail,  narrower  in  the  center  (about  2  earth  radii  thick)  and 
expanding  to  6  to  8  earth  radii  near  the  flanks  of  the  magnetosphere.  The 
change  in  thickness  of  the  plasma  sheet  can  be  explained  In  terms  of  charged 
particle  entry  and  magnetic  field  topology.  As  seen  in  Figures  3  through  5, 
substantial  amounts  of  plasma  enter  the  magnetosphere  only  within  about  20 
degrees  of  the  equator.  Quantitative  models  of  the  magnetospheric  magnetic 
field  show  that  field  lines  close  to  the  flanks  of  the  tail  are  more  dipolar 
in  shape  than  those  near  the  center  of  the  tail.  Thus  as  charged  particles  in 
the  plasma  sheet  cross  the  tail,  they  move  through  a  narrow  region  in  the 
center  of  the  tail  (see  Figure  7). 

Two  more  magnetospheric  plasmas  are  shown  in  Figure  7.  First  of  all  there  is 
a  boundary  layer  plasma  shown  just  inside  of  the  magnetopause  above  and  below 
the  lobe  regions.  This  boundary  layer  is  thought  of  as  an  extenstion  of  the 
magnetopause.  Technically,  however,  the  magnetopause  describes  that  region 
where  the  magnetic  field  abruptly  changes  magnitude.  The  boundary  l^yer 
plasma  can  be  explained  In  terms  of  our  entry  mechanism  as  consisting  of 
charged  particles  that  have  not  been  reflected  by  the  geomagnetic  field  but 
rather  have  gained  entrance  to  that  region  of  the  magnetosphere.  The  internal 
magnetic  field  topology,  however,  prohibits  these  particles  from  moving 
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further  into  the  magnetosphere.  Instead,  they  move  along  magnetic  field  lines 
down  the  tall  and  remain  confined  to  a  region  of  space  close  to  the 
magnetopause.  The  third  region  of  plasma  shown  In  Figure  7  Is  the  plasma 
mantle.  It  persists  just  Inside  of  the  boundary  plasma  and  at  mid  tall 
distances  has  a  thickness  of  about  1  earth  radius.  In  between  the  mantle 
plasma  and  the  plasma  sheet,  are  two  regions  referred  to  as  the  lobes  of  the 
tail  that  are  described  primarily  In  terms  of  the  magnetic  field  which 
persists  there.  The  lobe  regions  are  almost  devoid  of  plasma  because  they  are 
not  connected  magnetically  to  the  magnetopause.  That  is,  there  Is  no  way  for 
magnetosheath  particles,  even  once  they  have  gained  access  to  the 
magnetosphere,  to  move  into  the  lobe  regions.  In  the  terminology  of  cosmic 
ray  physics,  the  lobes  are  magnetically  forbidden  regions.  Another 
magnetically  forbidden  region  occurs  In  the  inner  magnetosphere  (see  Figure  6 
and  the  discussion  above). 

In  order  to  further  explain  the  existence  of  the  boundary  layer  and  the  plasma 
mantle,  it  Is  necessary  to  discuss  the  neutral  point  and  cusp  regions  of  the 
magnetosphere.  In  early  magnetospheric  magnetic  field  models,  it  was  found 
that  there  are  two  neutral  points  where  the  total  magnetic  field  vanishes. 

All  magnetic  field  lines  on  the  magnetopause  emanate  from  these  two  neutral 
points.  Later  it  was  observed  that  the  real  magnetosphere  does  not  possess 
neutral  points,  but  in  their  place  the  "day  side  cusps".  This  cusped  geometry 
is  shown  in  Figure  8  in  a  noon-midnight  meridional  cross  section.  Along  the 
magnetopause,  it  is  seen  that  entering  particles  are  free  to  move  over  a 
confined  region  that  constitutes  the  boundary  layer.  It  is  noted  that 
particles  which  enter  the  magnetosphere  deep  In  the  cusp  region  form  an 
expanded  plasma  layer  as  they  propagate  back  to  the  tall.  This  Is  the  plasma 
mantle.  The  reason  that  the  plasma  mantle  expands  to  about  one  earth  radius 
is  explained  simply  in  terms  of  the  change  In  the  strength  of  the  magnetic 
field.  In  order  for  the  charged  particles  to  maintain  their  magnetic  moments, 
their  gyroradll  must  increase  Inversely  proportional  to  the  strength  of  the 
magnetic  field.  Since  they  move  from  a  magnetic  field,  which  has  a  strength 
of  50  to  75  gammas  to  a  region  In  the  tail  where  the  field  has  weakened  to 
only  a  few  gamma,  the  gryo  radii  of  the  particles  Increase  from  about  100 
kilometers  to  a  substantial  fraction  of  an  earth  radius. 


It  Is  noted  that  all  of  the  work  reported  here  has  been  done  for  static 
magnetospheric  conditions.  We  realize  that  many  Interesting  magnetospherl c 
phenomena  require  an  understanding  of  the  magnetosphere's  dynamic  response  to 
changes  In  solar  wind  and  interplanetary  parameters.  In  fact  our  long  term 
goal,  stated  elsewhere  In  this  report.  Is  to  understand  and  quantitatively 
describe  the  time  varying  electric  and  magnetic  fields  which  reside  in  the 
magnetosphere  so  that  they  m^y  be  used  to  study  dynamic  plasma  behavior  in  the 
earth's  magnetospheric  "laboratory".  This  work  on  plasma  entry  should  help  us 
toward  that  goal.  It  is  interesting  to  note  that  particles  which  constitute 
the  boundary  layer  and  plasma  mantle  stay  in  the  magnetosphere  only  several 
minutes  since  they  flow  along  magnetic  field  lines.  This  is  in  striking 
contrast  with  the  particles  that  move  across  the  plasma  sheet.  Typically  it 
will  take  a  one  kilovolt  proton  with  a  small  pitch  angle  approximately  3  days 
to  move  across  the  tail.  This  difference  in  "magnetospheric  residence"  times 
has  important  implications  for  magnetosphere  dynamics  since  magnetospheric 
substorms  occur  typically  every  few  hours.  Thus  a  particle  moving  across  the 
plasma  sheet  will  probably  be  energized  and  deenergized  over  a  dozen  times 
during  its  journey  through  the  magnetosphere  or  it  may  be  injected  into  the 
ring  current  region  and  finally  precipitate  into  the  atmosphere. 

3.0  MAGNETOPAUSE  MICROPHYSICS 

In  a  sense  the  work  described  in  the  preceding  section  has  been  an  "aside". 

It  is  peripheral  to  our  primary  goal  and  yet,  like  our  earlier  excursion  into 
magnetic  fields,  we  feel  that  the  work  on  charged  particle  entry  into  a 
magnetically  closed  magnetosphere  is  an  important  subject  in  magnetospheric 
physics  in  its  own  right.  However,  we  have  continued  work  on  magnetopause 
microphysics  in  an  attempt  to  reach  our  goal  of  quanitlati vely  understanding 
the  formation  of  charge  separation  in  the  tail  region  that  supports  the  cross 
tail  electric  field.  We  have  found  this  to  be  an  exceedingly  complex  problem 
and,  in  fact,  one  unlike  any  others  we  have  worked  on  to  date.  This  is 
because  it  is  both  local  and  global  as  explained  below.  We  have  shown  above 
that  a  fraction  of  the  incident  protons  enter  the  magnetopause,  a  portion  of 
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them  moving  across  the  tail  to  form  the  plasma  sheet.  However,  the  protons 
enter  only  on  the  dawn  side  of  the  magnetosphere.  None  enters  on  the  dusk 
side  because  of  the  field  topology.  The  opposite  is  true  of  electrons;  they 
enter  on  the  dusk  magnetopause,  but  not  on  the  dawn  side.  If  all  things  were 
equal,  eventually  the  electrons  and  protons  woula  drift  across  the 
magnetosphere  forming  both  the  plasma  sheet  and  the  cross  tail  current  system 
that  supports  the  lobe  magnetic  field.  However,  there  remains  an  exciting 
problem  since  the  entry  fractions  for  the  protons  and  electrons  are  not  the 
same.  Thus,  even  for  a  charge  neutral  isotropic  plasma  existing  in  the 
magnetosheath  region,  more  protons  will  enter  on  the  dawn  side  than  electrons 
on  the  dusk  side.  We  conclude  therefore  that  the  work  described  in  the 
preceding  section  is  approximate  and  that  in  order  to  properly  understand  (1) 
charged  particle  entry  into  the  magnetosphere,  (2)  the  microphysics  of  the 
magnetopause,  and  (3)  the  formation  of  the  cross  tail  electric  field,  it  is 
necessary  to  develop  a  self  consistent  description  of  the  magnetopause 
magnetic  field  and  the  magnetic  field  Interaction  with  the  magnetosheath 
plasma.  The  solution  to  this  complicated  problem  would  provide  quantitative 
information  in  each  of  these  areas.  The  solution,  however,  is  difficult 
because  in  addition  to  requiring  the  self  consistent  description  of  the 
interaction  between  the  plasma  and  the  field  at  a  point  on  the  magnetopause, 
it  must  take  into  account  the  effects  of  the  electric  field  caused  by  charge 
separation  across  the  magnetosphere  due  to  differences  in  the  entry 
efficiencies  of  electrons  and  protons. 

In  an  attempt  to  understand  the  entry  problem  on  a  more  physical  basis,  two 
one-dimensional  models  of  the  magnetopause  were  developed.  A  description  of 
the  mathematical  aspects  of  the  models  is  given  in  Appendix  A.  Note  that  the 
narrative  and  equations  in  the  Appendix  were  done  on  one  of  our  word 
processors  that  we  are  just  now  learning  how  to  use.  It  provides  text  book 
quality  copy  (the  work  really  was  done  by  our  group). 

In  the  first  one-dimensional  problem,  an  MHD  macrophysical  model  was  developed 
which  led  to  the  classical  pressure  balance  equation.  This  model  was  used  to 
show  that  any  electrical  field  present  in  the  magnetosphere  must  be  very  large 


before  it  will  impact  pressure  balance  (on  the  order  of  1  volt  per  meter). 
Electric  fields  observed  in  the  magnetosphere  tend  to  be  in  the  millivolt  per 
meter  range  so  it  is  appropriate  in  the  pressure  balance  formalism  to  neglect 
the  pressure  caused  by  magnetospheric  electric  fields.  In  the  second 
one-dimensional  problem,  kinetic  theory  was  used.  Of  the  possible  solutions 
that  could  be  examined  only  one  was  carried  through  in  closed  form.  It  allows 
the  presence  of  the  electric  field  in  the  mangetosphere  but  does  not  permit 
particle  entry.  This  work  points  out  the  difficulty  of  determining  with 
mathematical  rigor  the  microphysics  of  magnetopause  structure.  Clearly  more 
work  must  be  performed  on  this  problem. 

4.0  PARTICIPATION  IN  THE  COORDINATED  DATA  ANALYSIS  WORKSHOPS 

Concurrent  with  our  woiic  on  magnetospheric  electric  fields,  we  have 
participated  actively  in  the  Coordinated  Data  Analysis  Workshops  (CDAW)  1,  2, 
and  6,  which  have  focused  on  the  dynamic  behavior  of  the  magnetosphere  by 
examining  specific  events.  The  following  paragraphs  contain  exerpts  of  work 
on  CDAW  described  In  our  December  81  proposal  to  your  office.  Since  the 
proposal  was  submitted,  the  first  of  the  CDAW  6  meetings  was  held  at  Goddard. 
We  also  describe  our  participation  in  that  meeting  as  it  relates  to  our 
studies  on  particle  entry  and  magnetospheric  electric  fields. 

While  we  were  working  on  the  entry  problem,  CDAW-1  and  CDAW-2  were  held,  both 
concerned  with  the  description  and  understanding  of  the  July  29,  1977  event 
while  at  CDAW-2  we  developed  the  first  quantitative  model  of  a  particular 
event.  The  magnetopause  current  strengths  and  location  were  determined 
directly  from  the  solar  wind  pressure.  Figures  9  and  10  show  the  solar  wind 
density  and  velocity  during  the  July  29  event.  The  pressure  balance  formalism 
developed  almost  20  years  ago  can  be  used  directly  to  calculate  the  size  of 
the  magnetosphere.  Figure  11  shows  the  standoff  distance  for  the  event.  Note 
that  early  in  the  event  the  magnetosphere  is  compressed  inside  of 
geosynchronous  orbit  and  that  later  the  magnetopause  relaxes  out  to  14  earth 
radii.  The  variation  in  the  strength  of  the  magnetopause  currents  through  the 
event  is  shown  In  Figure  12.  Thus  the  field  In  the  noon  quadrant  is 
accurately  modeled.  However,  in  the  antisolar  region  where  the  tail  currents 
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Figure  10 


MAGNETOPAUSE  STANDOFF  DISTANCE 


Figure  11 


STRENGTH  OF  MRGNETOPfiUSE  CURRENTS 


Figure  12 


become  important  the  agreement  between  observation  and  models  Is  not 
adequate.  Large  differences  between  observation  and  model  are  present  for 
both  events  In  this  region.  Thus  in  order  to  model  the  magnetospheric 
magnetic  field  in  the  most  dynamic  region,  the  tail,  more  work  Is  required. 
This  became  even  more  apparent  in  CDAW-6. 

At  the  CDAW-6  meeting,  another  pair  of  events  were  analyzed.  The  agreement 
between  the  model  and  the  magnetic  field  observed  by  the  SCATHA  satellite  is 
shown  in  Figure  13  and  seen  to  be  almost  amazingly  good  through  the  first  part 
of  the  event.  This  is  while  the  satellite  is  on  the  front  side  of  the 
magnetosphere.  Later  in  the  event,  however,  it  is  seen  that  large 
descrepancies  between  the  model  and  observations  appear.  In  an  attempt  to 
understand  the  source  of  the  discrepancies,  we  have  subtracted  the  model  field 
from  the  field  observed  by  the  GOES  2  and  3  satellites  at  synchronous  orbit 
near  midnight  (see  Figures  14  and  15).  This  residual  field  Is  that  portion  of 
the  observed  magnetic  field  that  cannot  currently  be  accounted  for  by  the 
dynamic  magnetic  field  models  we  have  used  for  the  CDAWs.  It  is  apparent  that 
this  discrepancy  cannot  be  explained  in  terms  of  a  ring  current  because  that 
would  also  impact  the  results  on  the  day  side  magnetosphere  where  agreement 
between  model  and  observations  is  excellent.  Also  the  magnetopause  currents 
are  too  distant  from  the  geosynchronous  satellite  locations  to  explain  the 
large  changes  in  magnetic  field  topology  observed  there.  Thus,  the 
discrepancy  is  caused  either  by  an  earthward  motion  of  the  plasma  sheet  and 
neutral  sheet  currents  (the  earthward  edge  of  this  current  system  would  have 
to  be  very  thin— a  small  fraction  of  an  earth  radius).  The  other  candidate 
current  system  is  one  that  is  linked  to  the  field  aligned  currents  which  close 
in  the  equatorial  regions  of  the  nightside  magnetosphere. 

Our  work  with  the  CDAW  is  thus  providing  us  more  information  on  the  location 
of  important  magnetospheric  current  systems  and  associated  plasma  phenomena  in 
electric  fields.  We  are  now  exploring  the  differences  between  the  observed 
night  time  disturbed  field  and  our  model  magnetic  field  and  attempting  to 
relate  them  to  our  work  on  the  physics  of  the  magnetopause  and  particle  entry 
into  the  magnetosphere. 
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GOES  -  3 
March  22,  1979 

Measured  field  ~  model  field 
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Figure  15. 


5.0  DISCUSSION  AND  CONCLUSIONS 

During  the  past  year,  as  described  above,  we  have  worked  on  the  microphysics 
of  magnetopause  structure,  the  entry  of  low  energy  charged  particles  Into  a 
magnetically  closed  magnetosphere,  and  have  been  active  participants  In  the 
6th  Coordinated  Data  Analysis  Workshop.  The  magnetopause  microphysics  problem 
is  very  interesting  In  Its  complexity.  In  fact,  it  should  be  of  general 
interest  to  the  plasma  physics  and  MHD  communities  since  it  is  an  example  of  a 
problem  that  requires  both  local  and  global  solutions  simultaneously.  This  Is 
because  the  structure  of  the  boundary  on  one  side  of  the  magnetosphere  Is 
dependent  on  charge  separation  and  resultant  electric  fields  caused  in  part  by 
events  on  the  other  side  of  the  magnetosphere.  It  Is  apparent  from 
discussions  with  investigators  currently  engaged  in  the  application  of 
3-Dimensional  magnetohydrodynamic  codes  to  the  magnetosphere  problem  that  the 
magnetopause  microphysics  problem  cannot  be  solved  with  an  MHD  formalism 
(because  of  the  problem's  local /global  nature).  We  have  continued  to  pursue 
the  problem  by  having  discussions  with  other  members  of  the  theoretical 
community  and  hope  to  understand  at  least  how  to  formulate  the  problem  in  the 
next  months. 

We  are  very  excited  about  our  particle  entry  work.  It  had  been  assumed  for 
many  years  that  only  very  energetic  particles  (solar  and  galactic  cosmic  rays) 
could  enter  Into  a  magnetically  closed  magnetosphere.  However,  It  Is  quite 
apparent  that  even  some  of  the  lowest  energy  particles  can  penetrate  the 
magnetopause  and  move  through  a  magnetically  closed  magnetosphere  owing  to  the 
finite  gradients  in  the  magnetospheric  magnetic  field  at  the  magnetopause. 

Most  of  the  entering  particles  are  constrained  because  of  magnetic  field 
topology  to  flow  along  magnetic  field  lines  in  the  anti  solar  direction. 
However,  near  the  center  of  the  magnetotail,  entering  charged  particles  can 
actually  move  across  the  tall  of  the  magnetosphere.  Thus  we  have  stated  that 
we  can  qualitatively  explain  the  existence  of  the  plasma  sheet,  plasma  mantle 
and  boundary  plasmas  in  the  magnetosphere.  To  date  we  have  not  performed 
detailed  calculations  of  plasma  densities  In  the  plasma  sheet,  etc.,  because 
of  the  lack  of  data  on  the  particle  distribution  function  for  the 
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magnetosheath  plasma.  Also,  our  work  on  particle  entry  has  been  restricted  to 
a  quiet  or  static  magnetosphere.  We  hope  In  the  near  future  to  extend  this 
work,  using  some  of  the  results  from  the  CDAWs,  to  Include  dynamic  processes 
occurring  in  the  magnetosphere  as  solar  wind  and  interplanetary  parameters  and 
thus  conditions  for  plasma  entry  change  during  "magnetospherlc  events". 

Our  work  at  CDAW-6  has  led  us  to  conclude  that  we  understand  the  day  side 
magnetosphere  very  well.  However,  there  muy  be  an  entire  current  system  that 
we  have  not  yet  Isolated  that  exists  In  the  near  tail  region  In  the 
magnetosphere  during  disturbed  times.  Most  likely  the  magnetic  residuals 
(differences  between  the  observed  and  model  values)  can  be  explained  in  terms 
of  the  extension  into  space  of  the  field  aligned  currents  which  flow  in  the 
ionosphere  at  auroral  altitudes.  We  plan  in  the  coming  months  to  develop  more 
realistic  input  parameters  to  our  existing  magnetic  field  models.  These  input 
parameters  will,  we  hope,  provide  information  on  the  dynamics  of  the 
magnetopause,  ring  and  tail  current  systems  that  are  based  on  a  physical 
understanding  of  particle  entry  into  the  magnetosphere.  The  ring  current  Is 
especially  interesting  since  particles  are  injected  sporadically  to  supply 
it.  Also,  particle  loss  rates  are  dependent  on  the  density  of  the  ring 
current.  Hopefully  these  phenomena  can  be  related  back  to  changes  in 
interplanetary  and  solar  wind  parameters.  We  expect,  however,  because  of 
energy  storage  within  the  magnetosphere,  that  for  both  the  ring  and  tail 
current  systems,  time  integrals  over  interplanetary  and  solar  wind  parameters 
will  have  to  be  developed. 
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I 


ONE  DIMENSIONAL  MODELS  OF  THE  MAGNETOPAUSE 


Introduction 

The  entry  of  charged  particles  through  the  magnetopause  is  ideally 
determined  when  these  particles  that  enter,  the  ones  reflected,  and  those 
drifting  through  from  within  the  magnetosphere  self-consistantiy  deter¬ 
mine  the  magnetopause  currents,  magnetic  and  electric  fields.  Classical 
magnetopause  theory  assumes  that  the  only  variation  of  the  various  field 
quantities  is  in  the  direction  perpendicular  to  the  magnetic  field  direction. 
This  assumption  will  yield  the  macroscopic  pressure  balance  equation  and 
specular  reflection  of  all  incident  particles.  The  presence  of  gradients  in 
the  magnetic  field  tangent  to  the  magnetopause  has  been  shown  to  allow 
particles  to  enter  a  magnetically  closed  magnetosphere  for  various  energies 
and  incident  directions.  This  calculation  was  not  done  self-consistently  and 
there  are  questions  what  form  the  pressure  balance  will  take  if  there  is  a 
net  entry  through  the  magnetopause.  Also,  if  electric  fields  are  formed  due 
to  the  differing  entry  efficiencies  of  protons  and  electrons,  how  do  these 
fields  effect  the  pressure  balance. 

We  assume  that  the  one  dimensional  geometry  used  in  classical  mag¬ 
netopause  physics  (see  Grad,  1961;  and  Hurley,  1963)  is  adequate  since  the 
gradients  of  field  quantities  parallel  to  the  magnetopause  are  much  smaller 
than  the  ones  perpendicular  to  the  magnetopause.  The  effect  of  these  small 
gradients  will  be  parameterized  by  the  current  perpendicular  to  the  mag¬ 
netopause.  Two  “solutions”  will  be  derived  here.  One  is  a  “macro”  solution 
which  is  the  equivalent  of  the  classical  pressure  balance  equation.  The  other 
is  a  “micro”  solution  which  illustrates  the  difficulty  of  doing  the  detailed 
kinetic  solution  even  in  the  simplest  cases. 

Assumptions 

We  make  the  one  dimensional  approximation  with  the  direction  per¬ 
pendicular  to  the  magnetopause  x,  and  the  B  field  in  the  z  direction.  The 
form  of  the  field  vectors  is 

B  =  B{x)i 
A  =  A(z)£ 

E  =  E(x)t 

where  we  have  the  relations  B  =  dAjdx  and  E  =  —d<p/dx  for  the  electro¬ 
static  potential  <p.  We  set  up  two  boundaries.  There  are  particles  incident 
from  x  =  —  oo  (i.e.  interplanetary  space)  with  prescribed  distribution 
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functions  Fi  defined  for  directions  with  vx  >  0  where  i  denotes  the  different 
species.  The  particles  incident  from  x  =  0  (i.e.  the  magnetosphere)  have 
the  prescribed  distribution  functions  G,  and  are  defined  for  directions  with 
vx  <  0.  Part  of  the  “micro”  solution  is  to  determine  the  total  distributions 
fi  for  —oo  <  x  <  0  that  are  due  to  the  boundary  fluxes  Ft  and  GL. 
We  assume  that  the  magnetic  field  is  zero  at  x  =  —  oo  and  increases  to 
some  unknown  value  at  x  =  0.  The  vector  potential  A  will  therefore  be 
monotonic  with  it’s  value  at  x  =  —  oo  set  to  zero.  The  potential  <p  is  also 
set  to  zero  at  x  =  —  oo.  We  will  assume  charge  neutrality  at  x  =  —  oo, 
and  if  possible,  at  x  =  0.  The  current  Jx  will  be  an  input  parameter.  A 
schematic  of  the  assumptions  is  shown  in  figure  1. 

The  equations  to  be  solved  are  the  Vlasov  equation  for  the  particles 
and  two  of  Maxwell’s  equations  for  the  fields 

v  •  dfi/d  r  (e»/m;)(E  +  vx  B)  ■  d/,/dv  =  0 

(d/dr)  •  E  =  —  d2<f>/dx2  =  (l/e0)  J  hdy  —  P/e o 

* 

f  ■  (d/dr)  x  B  =  —d2A/dx‘i  =  J  vyfxdv  =  VoJy. 

i 

In  the  remainder  of  this  section,  the  subscript  i  on  the  charge  e,  the  mass 
m,  and  the  distribution  functions  will  be  supressed  for  clarity  and  the 
summation  £  is  assumed  to  be  over  all  species  i.  The  distribution  function 
/  will  satisfy  Vlasov’s  equation  if  the  arguements  of  /  are  three  constants 
of  particle  motion.  We  use 

W  =$m2(v2x  +  v2y  +  v2z)  +  me(p 
P  —  mvy  +  eA 
Q  =  mvz 

where  W  is  particle  mass  times  the  total  energy,  and  P,Q  are  the  general¬ 
ized  momenta  in  the  y  and  2  directions.  We  will  transform  the  integrals 
over  velocity  space  into  integrals  over  {W,I*,Q)  space  according  (o  (lie 
relation 

,  dWdPdQ 


dv  — 


I  | 

d(  Va,V|,,Vs) 


where  the  Jacobian  is  given  by 


,  9{W,  P,Q) , 

9{vx,  Vy,  Vg) 


m4 \vx\ 


fi- 

t 


>  o 


The  particle  orbits  labeled  by  a  given  ( W,P,Q )  must  be  limited  ensure 
only  “allowed”  orbits  are  defined  by  the  distribution  function. 

Particle  Orbits 

Most  generally,  any  set  of  constants  (W,P,  Q)  at  a  position  z  v.  i?!i 
A(:r)  and  0(.r)  such  that 

2(W  —  me<p)  —  (P  —  eA)2  —  Q 2  —  rri2v2  >  0 

imply  that  the  particle  at  x  is  on  a  physically  valid  particle  orbit.  The 
eight  possible  forms  of  the  particle  orbits  is  shown  in  figure  2.  Note  that  no 
structure  in  the  y  or  z  directions  is  implied  in  this  figure.  For  all  particles 
incident  from  the  boundaries,  only  orbits  1-4  are  valid.  We  also  wish  for  the 
simplicity  of  having  all  orbits  connected,  so  oniy  cases  1-3  will  eventually 
be  allowed. 

To  avoid  cases  4  and  6-8;  we  must  not  allow  the  potential  defined 
by  tp  =  2 me<p  -}-  (P  —  eA)2  Q2  to  have  a  relative  maximum.  This  is 
the  condition  required  by  Grad  (1961)  derived  from  the  energy  equation, 
2 W  =  m2v 2  ip.  Setting  drp/dx  —  0  yields  P  =  eA  —  E/B.  Requiring 
that  d2tp/dx2  >  0  gives  us  the  condition  that 

(2 emE/B)dB/dx  +  2emp/eo  -f-  2e2B2  >  0, 

which  must  be  checked  after  a  kinetic  solution  is  derived.  Note  that  with  no 
electric  field,  the  result  reduces  to  Grad’s  restriction  that  A  is  monotonic; 
and  that  we  have  also  assumed  A  is  monotonic  so  that  B  can  only  go  to  zero 
asymptotically.  With  this  restriction  for  each  species,  the  possible  orbits 
are  reduced  to  the  4  shown  in  figure  3.  To  remove  case  5a,  we  require  that 
orbits  be  able  to  reach  either  the  left  or  right  boundaries. 

The  form  of  the  distribution  functions  limited  to  orbits  1-3  will  be 
derived  here.  Let  the  function  0(A,  <j>)  =  2 {W  —  em<p)  ~  (P  —  eA)2  —  Q2 
be  defined  where  the  (W,  P,  Q )  dependence  is  understood.  Let  /to  and  0O 
be  the  values  of  A  and  <p  at  x  =  0.  Then  the  distribution  funct  ion  /  is 
given  by 

f(W,P,Q-,A,(P)  =  F(W>P>Q)H[e(Q>0)]{2H[e(A)<l>)]  -  ff[©Mo,*o)]} 

+G(W,P,Q;/\o,0o)//r[e(/to,0o)]{2/f[e(A0)]  -  //[«((), 0)]} 

where  we  have  suppressed  the  dependence  of  /  on  (,40, 0())  and  II  is  lln 
Heaviside  step  function.  The  dependence  of  (7  on  (.40,  </>,,)  result'  from  < , 
given  in  terms  of  velocity.  Hecause  of  these  stop  functions,  the  iulcgration 
over  ( W,P,Q )  space  can  extend  over  the  entire  three  dimension  al  space 
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Physically  Passible  Orbits 

Figure  A-2 


Orbits  With  Restriction  On  * 


—  oo  <  W,  P,Q  <  oo  which  we  denote  by  5R3.  The  macroscopic  qualities 
number  density  n,  shielding  current  jy,  and  pressure  tensor  component  px  ,■ 
due  to  the  distribution  function  /  of  a  single  species  are 


dWdPdQ 


m3y/S(A,  (f>) 
e{P  —  eA)f(W,P,  Q',A,  <p) 
m*\J&(A,  <t>) 
v/e(A0) 


dWdPdQ 


mq 


f(W,  P,  Q;  A,  0) dWdPdQ. 


The  macroscopic  current  jx  of  a  single  species  is  given  by 


jx  ==  ~~~7  f  {F{W,  P,  Q)  -  G(W,  P,  Q;Aq,  0O)} 
ra4  Jx 3 

•  i/[©(0, 0)]tf[©(A0, 0o )]  dW  dP  d(J 


where  the  dependence  on  (Ao,0o)  for  this  as  well  as  the  other  macroscopic 
quantities  has  been  supressed. 

To  explain  the  use  of  the  step  functions,  note  that  m2v'j(A.(t>)  = 
©(A,0).  Using  the  x  —  — oo  boundary  flux  as  an  example;  the  particles 
that  come  from  this  boundary  must  have  an  initial  vx  >  0  so  ©(0,0)  > 
0.  If  they  reach  x  with  (A,  0)  =  (A{x),  (p(x)),  then  0(A,0)  >  0.  If  they 
reach  x  =  0,  they  will  not.  “reflect”  and  contribute  to  flux  with  v:r  <  0  r.o 
0(Ao ,0o)  >  0.  The  net  contribution  to  the  distribution  function  /  from 
the  left  boundary  flux  is  easily  seen  to  be 

F(M',/»,g)//[©(O,O)]{2//[©(A,0)]  -  //[©(AO,0„)]}. 


The  logic  for  determining  the  contribution  from  the  right  boundary  and  ju- 
is  similiar. 


Macro  Solution  -  The  Pressure  Balance 

Because  of  the  assumption  that  A  is  monotonic,  we  can  use  it  as  a 
coordinate  replacing  x.  Therefore,  the  electrostatic  potential  0  can  be 
written  as  <p  =  <p(A).  By  taking  the  total  derivative  of  the  pressure  pxx 
with  respect  to  A,  we  can  obtain 

dpxx{A,4(A))/tlA  —  jy  +  r.nP/lt. 

To  derive  this  result,  we  needed  to  use  several  of  the  relations  given  above 
Also,  in  taking  the  derivative  inside  the  integral,  we  get  terms  like  v/©<K©) 
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where  5  is  the  delta  function.  These  terms  integrate  out  to  zero.  Summing 
over  all  species  we  find 

p  =  £en 
Jy  —  52  Jy 
Jx  =  £  Jx 

Pxx  =  52  Pxx 

dPxx/dA  =  Jy-\-  pE/B. 

From  Maxwell’s  equation,  we  find  that 

—d2A/dx2  =  n0(dPxx/dA  —  pE/B ) 


which  we  integrate  from  A  =  0  (i.e.  x  =  —  oo)  to  A  (i.e.  2).  The  first 
term  integrates  to 


d2A/dx2(dA/dx)  dx 


—  f  %d(dA/dx)2 /dxdx  —  (dA/dx)2 /2  =  fP/2. 

Jx= - OO 


The  second  term  integrates  to 

rA 

M=0 


f  dPxx/dA  dA  =  Pxx  —  Pxx{  0, 0) 

•M=0 


where  /**((), 0)  is  the  value  at  x  ~  —  00.  The  last,  term  becomes 


[  pE/BdA  =  f  tpE/B)(dA/dx)dx 
•M=0  J  x=  —  oo' 

=  f  pEdx=  I  eoE(dE/dx)dx  =  eoE2/i 
J  x=  —  OO  J  x= — 00 


Adding  these  integrals  up  and  using  the  relation  po^o  =  1  /c2  where  c  is 
the  speed  of  light,  we  get  the  pressure  balance  condition 


Pxx{  x  =  —  00)  =  Pxx  -f  B2/  2p0  —  E2/2p0c2. 


The  last  term  gives  the  electric  field’s  effect  on  the  classical  pressure  balance. 
Unless  the  electric  field  is  very  large  (on  the  order  of  1  volt/meter  for 
the  earth’s  magnetopause),  the  electric  field  has  only  a  small  effect  on 
the  pressure  balance. 
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Micro  Solution  -  The  Boundary  Value  Problem 


The  solution  to  the  kinetic  problem  involves  two  of  Maxwell’s  equations 
which  are  coupled  through  the  current  Jy  and  charge  density  p.  At  (he 
boundary  x  =  0,  the  quantities  (A0, 0O)  are  unknown,  but  the  can  be  found 
from  the  known  conditions.  The  current  Jx  and  the  charge  neutrality  at 
x  =  —oo  are  assumed  known.  Since  these  quantities  are  functions  of 
(Ao,  4>o),  we  can  in  principle  solve  for  these  boundary  values.  The  equations 
to  be  solved  are 

—d2<j>/dx2  =  p(A,<l>)/e0 
—dPAjdx2  —  poJy(A,<p) 
where  the  boundary  conditions  at  x  =  —  oo  are 

A  =  0 
<t>  —  0 


and  the  boundary  conditions  at  x  =  0  (i.e.  (Ao,0o))  are  found  from 

P(  0,0)  =  0 
Jx  =  •/known 

which  are  two  equations  for  the  two  unknowns  (A|),0»). 

Micro  Solution  -  A  Simple  Example 


To  illustrate  the  above  micro  solution,  we  solve  a  simple  example  that 
allows  no  particle  entry,  but  does  form  an  electric  field.  We  assume  that  a 
beam  of  mono-energetic  particles  is  incident  normally  on  the  magnetopause. 
The  distribution  functions  for  incident  ions  and  electrons  are 


Fi  = 


mfy/2  rcpc 

2  y/WiPiQ* 


6(W/Wt-l)6(Q/Qt)6(P/Pt) 


and 


K 


m?\/2  Woo 
2,/W;.PrQ'. 


t>{w/w,  -  i VHQ/Qrwnfrr) 


where  W,,W,  >  0.  Letting  r.  be  the  positive  electron  charge,  the  integral 

Jx  is 


h 


— — {(  y/wif  m,i)H\2{Wi  -  em,0o)  - 

x/2 


-(\/vn/me)tf  [2(We  +  emt<p0)~  e2A^]}. 
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We  assume  Jx  is  zero  so  for  general  Wi  and  We,  we  have 


2 (W»  —  em,0o)  —  e2A2  <  0 

2 {We  +  eme<p 0)  —  e2A2  <  0. 
The  current  Jy  is  given  by 


(*) 


y/WiH[2{Wi  -  emit)  -  e2A2] 
rrii\/2{Wi  —  emit)  —  PA1 

^  y/WeH[2{We  +  eme<t>)  -  e2A2}  | 
me\/2{We  4-  erne<p)  —  e2A2  ) 


when  we  apply  the  above  condition  at  x  =  0.  In  fact,  we  can  let  the 
inequalities  in  (*)  be  equalities  so  that  the  current  Jy  will  exist  up  to  x  = 
0.  The  charge  density  is  given  by 


enc 


f  SW  H\2(Wj  -  emit)  -  e2A2] 
A  s/2{Wt  -  ermt)  -  e2A2 


y/2WeH[2{We  4-  eme<t>)  -  e2A2]  | 
y/2(We  4-  eme$)  —  e2A~  / 


where  the  step  functions  will  be  equal  to  1  in  the  range  —  oo  <  x  <  0  by 
the  equality  condition  in  (*). 

The  charge  density  derived  above  can  be  set  to  zero  so  that  a  quasi¬ 
neutral  relation  between  <f>  and  A  is  obtained.  The  result  is  that 

<P  =  MA/Ao )2 


where 


00 

Aq 


1  Wj-We 

e  mi4-me 

2  mjWe  4-  mc'.Vj 
e2  mi  4-  me 


are  the  solutions  to  the  equalities  in  (*).  The  only  differential  equation  left 
to  solve  is  Maxwell’s  equation 


d2A/dx2  =  —  HoJy  = 


/*oe2noo  f  +  _1_^ 
y/\  -  A2/A^Ue  rnj 


The  change  of  variable  from  x  to  a  —  A/Ao  will  transform  the  equation 
to 

iB2/ia  =  -  gjg= 

Vl-  a2 

where  B%  =  2Ale2n00(iQ(l/me  +  1  /m*).  The  solution  of  this  differential 

equation  is  _ 

B 2  =  B2(  1  —  Vi  —  a2) 

where  the  value  of  B2  at  x  =  0  (i.e.  a  =  1)  is  in  fact  B2  defined  above. 
Since  the  x  derivative  of  A  is  B,  the  solution  for  B 2  gives  us  another 
differential  equation  to  solve.  We  can  manipulate  it  into  the  form 


dx/da  — 


Ao 


Bq\Ji  —  y/T— 

which  has  the  normalized  solution 


f*. = 1) + v-2  u 


Figure  4  shows  the  dependence  of  normalized  A,  B,  and  <j>  on  normalized 
distance  for  this  simple  problem. 


DISCUSSION 

The  kinetic  example  above  was  the  simplest  magnetopause  problem  of 
interest  that  could  be  solved  in  closed  form.  In  this  quasi-neutral  solution, 
we  could  throw  out  one  differential  equation  and  luckily  the  other  could 
be  solved  analytically.  There  is  no  current  through  the  magnetopause  lml 
there  is  an  electric  Held  whose  sign  is  determined  by  the  particle  with 
the  larger  incident  momentum,  \Pi W.  Because  these  particles  are  mono 
energetic,  there  are  only  three  discrete  currents  Jx  possible.  A  more  realisl.it 
problem  with  incident  beams  having  a  continuous  energy  distribution  would 
allow  a  continuum  of  Jx  currents  to  be  specified;  however,  this  problem  will 
be  very  difficult  to  solve  if  it  is  possible  to  solve  at  all. 

The  kinetic,  micro/solution  approach  is  seen  to  be  difficult,  but  it  can 
give  us  all  the  details  of  the  solution  if  we  want  them.  If  all  we  want  is 
the  potential  <po  inside  the  magnetopause,  we  will  not  need  to  solve  the 
differential  equations.  The  incident  distribution  functions  must  still  be 
known  and  the  coupled  integrals  for  p(x  =  —  oo)  and  Jx  solved  for  (A0,  <2>o)- 
The  macro/solution  will  not  yield  the  electric  field  since  it  will  be  a  small 
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Figure  A-4 


residual  of  the  classical  pressure  balance.  Because  the  pressure  balance  is 
essentially  a  MHD  equation,  the  MHD  simulations  can  not  be  expected  to 
yield  an  electric  field  formed  by  particle  interaction  at  the  magnetopause 
even  if  the  resolution  near  the  magnetopause  was  adequate. 

A  possible  way  to  solve  the  magnetopause  problem  with  out  having  to 
solve  Maxwell’s  equations  would  be  to  use  B  field  gradients  to  determine 
current  Jx  and  solve  for  (Ao,  <p o)  as  described  above;  substitute  the>e  values 
in  the  integrals  for  Pxx  inside  and  outside  the  magnetopause;  and  use  the 
pressure  balance  to  determine  B  inside  the  magnetopause.  The  distribution 
functions  for  particles  incident  from  inside  the  magnetosphere  are  still 
needed.  These  nave  entered  from  another  position  on  the  magnetopause, 
so  a  global  solution  oe  the  state  of  the  entire  magnetosphere  is  necessary. 
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